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ERROR REDUCTION PROGRAM* 

A PROGRESS REPORT 

Saadat A, Syed 

Engineering Division; Pratt & Whitney Aircraft Group 


The combustion chamber of the aircraft gas turbine has to be designed to satisfy 
the stringent dual engine requirements of reduced fuel consumption and increased 
durability. Combustor development based on conventional methods to meet these 
requirements is proving to be increasingly costly and time consuming. A mathe- 
matical model of the combustcr would be very helpful In reducing cost and time of 
the design cycle as it would permit most of the work to be done on a computer. Such 
models are already being used in a limited manner in the industry. More widespread 
use of these models will be accelerated by the removal of some of the known 
deficiencies in the present codes. One of these deficiencies is the numerical 
error, or numerical diffusion, that can be introduced into the calculations under 
certain conditions. This numerical error is due to the finite differencing scheme, 
usually upwind or hybrid, used in these codes. 

The objective of the Error Reduction Program is to evaluate available finite 
difference schemes for minimum numerical diffusion, and to incorporate the best 
available scheme into a three-dimensional combustor performance code. 

As a first step toward this end, five schemes (refs. 1 through 5) were selected 
for initial evaluation. Some of these schemes were combined with bounding schemes 
to eliminate physically unrealistic undershoots and overshoots. Two basic criteria 
were used to judge these schemes: they should be conservative, and should produce 
solutions that exhibit no extraneous maxima or minima (boundedness property). The 
accuracy of the schemes was evaluated by performing the truncation error analysis, 
and running one- and two-dimensional test cases and comparing the calculated 
solutions against the exact solutions. Based on this evaluation, two schemes were 
selected: QUDS, Quadratic Upstream Differencing Scheme, and BSUDS2, Bounded Skew 
Upstream Differencing Scheme Two. The first scheme was proposed by Leonard (ref. 
1), and the second scheme is the one proposed by Raithby (ref. 2), which is bounded 
by a new bounding scheme. 

The selected two schemes were coded into a two-dimensional computer code, 2D-TEACH, 
and their accuracy and stability were evaluated by running several test cases (refs. 
6 through 8). It was found that BSUDS2 was more stable than QUDS. It was also 
found that the accuracy of both schemes is dependent on the angle that the 
streamlines make with the mesh, QUDS being more accurate at smaller angles and 
BSUDS2 being more accurate at larger angles. However, both schemes, at all angles, 
were more accurate than the existing hybrid scheme. BSUDS2 was selected to be 
extended into three dimensions, primarily because It was more stable. This scheme 
is currently being incorporated into a three- dimensional code, 3D-TEACH. 


* This work was conducted under NASA Contract NAS3-23686. 
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NEED FOR ERROR REDUCTION OBJECTIVE 



Conviction and diffusion of a scalar stop m constant proparty 
uniform wh toned flow Pa a Pa y 90 


To identify and incorporate best available 
finite difference error reduction scheme 
into a 3D combustor performance code. 
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QUAUTV DISCRETIZATION SCHEMES ; 

EVALUATED 


• Agarwal Differencing Scheme (ADS) 

• Cuadratic Upwind Differencing Scheme (QUDS) 

• Skew Upwind Differencing Scheme (SUDS) 

• Spline/Hermetian Schemes 

• Cubic Spline Scheme (CSS) 

• Glass and Rodi Hermetian Scheme (GRHS) 

• Flux Blending Schemes 

• Bounded-One 

• Bounded-Two 


SCHEMES SELECTED 


• QUDS — Quadratic Upwind Differencing Scheme 

• BSUDS2 — Bounded Skew Upwind Differencing 

Scheme Two 
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WHAT IS BOUNDEDNESS 

In the absence of sources the value at node P should 
be bounded by the values at surrounding nodes. 
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ORIGINAL PAGE 19 
OF POOR QUALITY 


QUADRATIC UPWIND 
DIFFERENCING SCHEME 
(QUDS) 

• Finite volume method 

• Upwind biased quadratic interpolation for convection terms 

• Central differencing for diffusion terms 

• Employs extended nine point modecule 

= ~ V*^VVW + + 3/ *^P u > 0 

= - + 3 /«^W' u < 0 

• Scheme is conservative 

• Coefficients can become negative 

• Solution can become unbounded 



i 


i 


i 

\ 

r 

i 

i 

i 


BOUNDED SKEW UPWIND 
DIFFERENCING SHCEME TWO 
(BSUDS2) 

• Finite volume method 

• Central differencing for diffusion terms at all Peclet 

numbers 

• Central differencing for convection terms for absolute 

Peclet numbers less than two 

• Blending of upwind and skewed , t 

upwind differencing for Peclet 
numbers greater than two 

• Employs compact nine-point 1 1 . \ 

molecule 'j- . . 

! 

• Scheme is conservative * -.j 1 

• Produces bounded solutions 
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ORIGINAL PAGE 19 
OF POOR QUALITY 


FIRST LAMINAR FLOW TEST CASE 

QUDS more accurate over all 
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SECOND LAMINAR FLOW TEST CASE 
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STAGNATION STREAMLINES FOR SECOND 
LAMINAR FLOW TEST CASE ON 40 x 20 MESH 



STAGNATION STREAMLINES FOR 
SECOND LAMINAR FLOW TEST CASE 
ON 78 x 40 MESH 
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FIRST TURBULENT TEST CASE 
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Inlet boundary conditions 


Mean velocity 

Temperature 

Pressure 

Boundary layer thickness 
Turbulence intensity 


60 1 1 ft/sec 
59'F 
14 7 psia 
0 33 inch 
0 003 


118 32 m/s) 

(41 2"C) 

(101 35kN/m 2 ) 
(8 4 mm) 
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CALCULATED REATTACHMENT POINTS 
(MEASURED 7 STEP HEIGHTS) 


Case 

HYBRID 

BSUDS2 

QUDS 

Dumb (15x17) 

2.4 

2.25 

2.28 

Coarse (26x29) 

5.2 

5.4 

5.5 

Fine (50x56) 

5.7 

5.9 

Unstable 

Fine adjusted 
(74x53) 

5.8 

5.8 

Unstable 


- 
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SECOND TURBULENT TEST CASE 



SECOND TURBULENT TEST CASE 
WITH SWIRL 


Inlet boundary conditions; 

* Measurements at x = 5mm were used 

• This plane was considered as inlet plane 


SECOND TURBULENT TEST CASE 
WITH SWIRL 
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Stations at which comparisons have been made 
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SECOND TURBULENT TEST CASE WITH SWIRL 
COMPARISON OF 
AXIAL VELOCITY PROFILES 

C oars# grid (40x35) 
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BSUDS2 
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No difference in accuracy for fine grid (47x43) 
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SECOND TURBULENT TEST CASE WITH SWIRL 
COMPARISON OF 
TANGENTIAL VELOCITY PROFILES 

Coarse grid (40x35) 

HYBRID 

QUDS 

BSUD82 
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QUD8 seems more accurate 
No difference In accuracy for fine grid (47x43) 
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SINGLE CELL CALCULATIONS 
OF CONVECTION OF NORMAL DISTRIBUTION 


^(xl - exp ( VjW A*)M 
I 7 ? x/A,A # - A/cr) 



0, degree 

Influence of the flow angle on the error 
for various difference schemes 
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OBSERVATIONS 

• QUDS unstable — needs a better solver 

• No significant differences between HYBRID, QUDS 

and BSUDS2 for intelligent fine grids 
(1000-1500 nodes) 

• QUDS more accurate than BSUDS2 most of the time 

• BSUDS2 more accurate for second laminar test case — 
very strong streamline curvature 

• Flows with swirl more sensitive to difference schemes 
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